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We study the production of J/ip + cc in e + e~ annihilation through two 
virtual photons. The cross section is estimated to be 23 fb at = 10.6 
GeV, which is smaller by a factor of six than the calculated cross section 
for the same process but through one virtual photon. As a result, while the 
annihilation into two photons may be important for certain exclusive pro- 
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duction processes, the big gap between the inclusive production cross section 
a(e + e~ — > J/ip + cc) ~ 0.9pb observed by Belle and the current nonrelativistic 
■ QCD prediction of ~0.15pb still remains very puzzling. We find, however, as 

£f~) • the center-of-mass energy increases (y/s > 20 GeV) the production through 

two virtual photons e + e~ — > 27* — ► J/ip + cc will prevail over that through 
one virtual photon, because in the former process the photon fragmentation 
into J/ip and into the charmed quark pair becomes more important at higher 



Q-i. energies. 



PACS number(s): 12.40.Nn, 13.85.Ni, 14.40. Gx 



Charmonium production is one of the important processes to test quantum chro- 
modynamics (QCD) both perturbatively and non-perturbatively. Because of the 
simpler parton structure involved, which will reduce the theoretical uncertainty, 
charmonium production in e + e~ annihilation is expected to be more decisive in 
clarifying the production mechanisms of heavy quarkonia. The two B factories with 
BaBar and Belle are collecting huge data samples of continuum e + e~ annihilation 
events, which will allow us to have a fine data analysis for charmonium production. 
Recently, Charmonium production in e + e~ annihilation has become more interesting 
and puzzling, because of the large gap between the Belle measurements)]] E] an d the 
theoretical calculations for both inclusice[3lll|31in] and exclusive [HE] charmonium 
production via double cc pairs based on nonrelativistic QCD (NRQCD). 

For the inclusive processes, Belle has reported a measurement on the J/ip produc- 
tion in e + e~ annihilation at \/s = 10.6 GeV [HE], and found that a very large fraction 
of the produced J/ip is due to the double cc production in e + e~ annihilation [T] 

a(e + e~ -> J/ipcc)/a(e + e- -> J/ipX) = 0.59±g;i| ± 0.12, (1) 
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which corresponds to |2| 

a(e + e~ -> J/^cc) » 0.9p6. (2) 

In contrast, the predicted values for the cross section in NRQCD (the color- 
octet contribution is negligible for this process) are much smaller than the data[3, 
HJIH]. In a recent analytical calculation and numerical estimation for the inclusive 
charmonia production including all S-wave, P-wave and D-wave states via double cc 
in NRQCD 0, we find 

a(e + e- -> 7* -> J/^cc) « 0.15p&. (3) 

This value is consistent with other previous calculations, including those obtained 
based on the quark-hadron duality hypothesis [9 , but smaller than the Belle data 
by about a factor of six 1 . 

For the exclusive processes, the Belle measurement [T] for the cross section of 
e + e~ — > J/ip + rjc is about an order of magnitude larger than the NRQCD calculation 
for e + e~ — > 7* — > J/ip + r/cjH IE]- In order to solve the problem, calculations for the 
exclusive double-charmonium production from e + e~ annihilation into two virtual 
photons are performed |in|, and it is pointed out that the cross section for e + e" — > 
27* — > J/ijj + J/ij) is larger than that for e + e~ — > 7* — > J/ip + f] c by about a factor of 
3.7, despite of possible uncertainties due to the choice of input parameters [TT]. This 
is a interesting result since it indicates that the e + e~ annihilation through two virtual 
photon fragmentation may make important contributions to certain processes, and 
it might substantially reduce the discrepancy between experiment and theory for the 
exclusive process e + e~ — > J/ip + rj c (note that it is essential to check experimentally 
whether e + e~ — ► J/ip + J/ip is largely misidentified with e + e~ — > J/ip + rj c ). 

In this situation, it is necessary to examine the contribution of e + e~ annihilation 
through two virtual photons to the inclusive production of J/ip. In the following we 
will calculate the complete 0(a 4 ) color-singlet inclusive cross section for e + e" — > 
27* — > J/ipcc, and compare the production rate through two virtual photons with 
that through one virtual photon, to see whether the annihilation through two virtual 
photons can decrease the discrepancy between the Belle measurement on J/i/jcc 
production and the calculations based on NRQCD. 

Following the NRQCD factorization formalism, color-singlet scattering ampli- 
tude of the process e _ (px) + e + (p 2 ) —* 7* — > cc( 25+1 -Lj)(p) + c(p 3 ) + c(p 4 ) in Fig. [T] 
is given by 

A{e-{ Pl ) + e + (p 2 ) - cc( 2S+1 Lj)(p) + c(p 3 ) + c(p 4 )) = ^EEE 

L Z S Z jk 

x(si,s 2 I SS Z )(LL Z ;SS Z I JJ z )(3j;3k \ 1) 

xA(e-( Pl ) + e + (p 2 ) - Ci (|; Sl ) + c fc (|; s 2 ) + q(|; s 3 ) + c,(|; s 4 ))(L = S). (4) 

where cc( 2 +1 Lj) is the intermediate cc pair which is produced at short distances and 
then evolves into a specific charmonium at long distances. (3j; 3k | 1) = 5jk/ '\fN~c , 

x The numerical value obtained for this process in 3 should be multiplied by a factor of 3. 
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Figure 1: Feynman diagrams for e + e — > 27* —>■ J ftp + cc. 

(si;s2 I S'S'z) , (LL Z ;SS Z \ JJ Z ) are respectively the color-SU(3), spin-SU(2), and 
angular momentum Clebsch- Gordon coefficients for QQ pairs projecting out appro- 
priate bound states. For J/ip the coefficient Cl can be related to the radial wave 
function of the bound state and reads 

Cs = i I Rs(0) | 2 • (5) 
The spin projection operator can be defined as [T2*] 

PssAPW) = ^(si,S2\SS z )v(~ + q;s 1 )u(~-q;s 2 ). (6) 
We write the spin projection operator which will be used in the calculation as 

Pis z (P,0) = ^AS z )(fi + M), (7) 

where M is the mass of the charmonium, which equals to 2m c in the nonrelativistic 
approximation. 

The amplitude for the upper Feynman graph in Fig. Q can be written as 

Tr[P 15z (p, 0) 7a ]lZ(p 3 )7 M t;(p4) a , * , 2 , (8) 

3 



100 n 



80- 



q 60- 

w 
o 
o 

^ 40- 

D 
"O 



20- 



-1.0 -0.5 0.0 0.5 1.0 

cose 



Figure 2: Differential cross sections of e + e — ► 27* — > J/ip + cc (solid line) and 
e + e~ — > 7* — > J/-0 + cc (dotted line) as functions of the scattering angle of J/^>. 




and the amplitude for the corresponding flipped graph is denoted as M%. The 
amplitude for the lower Feynman graph in Fig. ^ reads 

AjT 4 X X^/TT CC I TT\~f \ fa+ £/ 2 ~ fa + m e / x 

M3 = Vl^ S ( ' JJ>(P2)7 ^ + P/2 - V,Y - mf aM(Pl) 



1 

(p 4 + p/2) 2 (P3+p/2) 2 ' 



m(P3)7/.-Pls,(p, 0)70^^4) ^ r^TT^TI 1 „/o\a > ( 9 ) 



and the amplitude for the flipped graph is denoted as M 4 . 

The calculation of cross section for e~ + e + — > 27* — > J/^ + cc is straightforward. 
The differential cross section can be written in the form 

da{e + + e -> 2 7 * -> + cc) = 3C s a 4 ^ 2 

where cZQ represents the elements of the quadruple integral related to the final state 
phase space (see, e.g. Ref. [3]), and | M | 2 = jJ2 P oi,coior I M± + M 2 + M 3 + M 4 
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is 



the unpolarized module squared of the amplitude. For simplicity we will not write 
down their lengthy expressions here. 

With Eq. (JTUJl we can evaluate the inclusive cross sections for J/i/j production 
from e + e~ through two virtual photons. The input parameters used in the numerical 
calculations are the same as in ref. (i 

m e = 0, m c = l.5GeV, a = 1/137, (11) 

I r s (q) | 2 = 0.81GeV 3 |H|. (12) 
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Figure 3: Ratio of the fragmentation contribution to the total cross section as the 
function of the center-of-mass energy. 



Now we give the numerical result at the Belle energy yfs = 10.6 GeV: 

a(e + + e~ -> 2 7 * -> J/ip + cc) = 23fb. (13) 

It is well known that e + e~ — > 27* — > J/ip + cc is a pure electromagnetic process 
(except for the hadronization of the quark pair at long distances), while e + e~ — > 
7* — > J/ip + cc involves both electromagnetic and strong interactions. For a naive 
order of magnitude estimate, the ratio of the production rate of the former to the 
latter would be proportional to a 2 /a 2 , but the photon fragmentation into J/ip and 
into the charmed quark pair can substantially enhance the former. To see the 
role of photon fragmentation, we plot the differential cross sections of e + e~ — > 
27* — > Jjijj + cc and e + e~ ^7* — > J/ijj + cc as functions of the scattering angle 
of J/ip at = 10.6GeV in Fig. |2j (Here for the one-photon process we choose 
a s =0.26.) One can see that the small angle J/ip production in e + e _ — > 27* — > 
J/i/j+cc is dominant. This indicates that most of the J/ijj production comes from the 
photon fragmentation (corresponding to the upper graph in Fig. ^ with one photon 
fragmenting to J/ip and the other fragmenting to a charm quark pair). Indeed, 
our calculation shows that at y/s = 10.6GeV, the contribution from fragmentation 
graphes is about seventy-two percent of the total cross section. In Fig. El we show 
the ratio of the fragmentation contribution to the total cross section as a function 
of the e + e~ center-of-mass energy in e + e~ — > 27* —>■ J/tp + cc. It is clear that 
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Figure 4: Cross sections of e + e~ — > 27* — > J/ijj + cc (solid line) and e + e~ — > 7* — > 
J/ip + cc (dotted line) as functions of the center-of-mass energy 

the photon fragmentation becomes more and more dominant as the center-of-mass 
energy increases. This is in agreement with the observation in the J/ip + J/ip 
exclusive production through two virtual photons [TP]. 

From the above discussions we have seen the importance of the photon frag- 
mentation to J/ip as well as to the charm quark pair in the two-photon process 
e + e~ — > 27* — > J/ip + cc. A even more crucial result is that at high e + e~ center- 
of-mass energies the contribution through two virtual photons will prevail over that 
through one virtual photon in the production of J/ipcc. The reason lies simply in the 
fact that the virtuality of the photon in the two-photon process can be as small as 
4m 2 , whereas it is as large as s, the center-of-mass energy squared in the one-photon 
process. In Fig. 0]we show the cross sections of e + e~ — > 27* — > J/ip + cc (solid line) 
and e + e" — > 7* — >• J/ip + cc (dotted line) as functions of the center-of-mass energy 
y/s. We see clearly that the cross section for e + e~ — > 27* — > J/ip + cc decreases very 
slowly, whereas that for e + e~ — > 7* — > J/ip + cc decreases rapidly as yfs increases. 
At y/s = 20 GeV, the two photon process becomes to prevail over the one photon 
process. 

However, unfortunately, at the Belle energy y^s = 10.6 GeV, since the enhance- 
ment effect due to the factor s/m 2 c is not large enough as compared with the sup- 
pression factor a 2 /a 2 s1 we find a(e + e~ — > 27* — > J/ip + cc) = 23 fb, which is still 
much smaller than cx(e + e~ — > 7* — > J/ip + cc) = 148 fb|6J, and therefore is negligible. 




6 



In summary, we have calculated the complete O (a 4 ) color-singlet inclusive cross 
sections for J / ipcc production from e + e~ annihilation into two photons. Due to the 
suppression factor of a 2 /a 2 s) at the e + e~ center-of-mass energy s/s = 10.6 GeV the 
cross section of this process is smaller by about a factor of six than that from e + e~ 
annihilation into one photon. We then conclude that while the e + e~ annihilation into 
two photons could be helpful in solving the puzzle for the exclusive J / iprj c production, 
it can do very little to reduce the big gap between the observed inclusive production 
cross section of cr(e + e~ — > J/ip + cc) ~ 0.9 pb and the current NRQCD predictions 
of about 0.15 pb. This puzzle still needs to be explained with new theoretical 
considerations. We find, however, as the center-of-mass energy increases (-^/s >20 
GeV) the production through two photons e + e~ — > 27* — > J/ip + cc will prevail over 
that through one photon, because in the former case the photon fragmentation into 
J/ip and into the charmed quark pair becomes more important at higher energies. 
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